Replenish and relax: explaining logarithmic annealing in disordered materials 
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Fatigue and aging of materials are, in large part, determined by the evolution of the atomic-scale 
structure in response to strains and perturbations. This coupling between microscopic structure 
and long time scales remains one of the main challenges in materials study. Focusing on a model 
system, ion-damaged crystalline silicon, we combine nanocalorimetric experiments with an off-lattice 
kinetic Monte Carlo simulation to identify the atomistic mechanisms responsible for the structural 
relaxation over long time scales. We relate the logarithmic relaxation, observed in a number of 
systems, with heat-release measurements. The microscopic mechanism associated with logarithmic 
relaxation can be described as a two-step replenish and relax process. As the system relaxes, it 
reaches deeper energy states with logarithmically growing barriers that need to be unlocked to 
replenish the heat-releasing events leading to lower energy configurations. 
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The importance of defects in determining material 
properties and evolution can hardly be overstated, and 
yet understanding their evolution over experimental 
timescales remains a challenge. This problem affects ma- 
terials ranging from somewhat simple crystalline systems 
to glasses, alloys and cements and is at the roots of phe- 
nomena affecting materials reliability, aging and fatigue. 
Over the years, their evolution has been characterized 
by various models based on continuous, mesoscopic or 
atomic scale models and generally fitted to experimental 
data [2 [2]. Direct comparison between models, experi- 
ments, and atomic-scale simulations, however, have been 
limited because of the often extensive time scale differ- 
ence between the three, leaving considerable space for 
interpretation in the fitting procedure even when describ- 
ing apparently simple relaxation processes such as loga- 
rithmic relaxation observed in polymer glasses, for exam- 
ple [3l|4], and in ion-implanted amorphous and crystalline 
systems [5l[6]. Here, we select ion implantation induced 
disorder in monocrystalline silicon (c-Si), a controlled ap- 
proach for perturbing samples that can be closely repro- 
duced numerically [7, 8 . Combined with calorimetric 
measurements, ion implantation is an ideal technique for 
studying the kinetics and thermodynamics of complex 
structures evolution in a systematic way. 

Many models have been proposed for describing the 
accumulation and annealing of ion-induced disorder in 
materials. The Frenkel pair model, for example, de- 
scribes the disorder in terms of isolated vacancies (V) 
and interstitials (I) [9 ; in contrast, MD simulations 
show that the majority of defects induced by the im- 
pact of a single ion of a few keV are found in the form 
of heavily-damaged zones commonly called amorphous 
pockets (APs) [7, 10, llj that anneal in steps, with a wide 
range of activation energies that depend on the details of 
each APs' interface with the crystal [T", as observed ex- 



perimentally [12]. A popular model to account for dam- 
age accumulation proposed by Marques and coworkers 
[13j describes APs as an assembly of hound defects corre- 
sponding to IV-pairs [M]. While never detected directly 
in experiments, this mechanism was observed in various 
forms in short-time simulations |6l |T5l |T6] . As shown be- 
low, this model does not describe nanocalorimetriy (NC) 
measurements correctly and we are still seeking the cor- 
rect atomic-scale relaxation model for implanted semi- 
conductors. 

In this Letter, we develop such a model by coupling 
the NC measurements of the annealing of low-energy, 
low-fluence implantation-induced disorder in c-Si with 
results from simulations performed using the recently 
proposed kinetic Activation- Relaxation Technique. K- 
ART is an off-lattice, self-learning kinetic Monte Carlo 
algorithm [I8l [19] that makes it possible to follow the 
time evolution of large implanted simulation boxes over 
a timescale of up to a second or more, taking into account 
the full complexity of activated events and their associ- 
ated long-range elastic effects over timescales comparable 
the experiments. Results show that the broad, feature- 
less heat release as a function of temperature measured 
by NC over several hundreds of Kelvin is well reproduced 
by the k- ART simulations. Analysis of the atomistic sim- 
ulations allows us to understand the microscopic origin 
of this logarithmic relaxation. 

NC measures heat release or absorption as a function 
of temperature in thermodynamic and kinetic processes 
occurring at the nanoscale. The technique, described in 
Refs. [20l|2T], has been used to investigate phenomena 
ranging from melting point depression in nanoparticles 
[22] to disorder annealing in amorphous Si (a-Si) [5 and 
polycrystalline Si (poly-Si) [23 . The device consists of a 
low-stress, 250 nm thick SiN^^ membrane supported by a 
300 jivd thick silicon frame. On top of the membrane, a 



Pt strip is used to measure the temperature and heat the 
sample. The new device used here [24 features, under- 
neath the membrane, a 330 nm thick c-Si strip, ahgned 
with the heating strip. Heating rates reaching 10^ K/s 
convert smah amounts of heat released during disorder 
annealing into measurable power. 

For these experiments, the sample and reference NCs 
were in contact with a sample holder either maintained 
at room temperature (RT) or cooled with liquid nitrogen 
(LN). In the latter case, the membrane, surrounded by 
a thermal shield, reached 110 K. The silicon layer was 
implanted with Si~ at a fluence of 0.02-0.1 Si/nm^ and 
an energy ranging from 10 to 100 keV. The current was 
measured and integrated with systematic error <20%, al- 
though the relative uncertainty between different exper- 
iments within a series of implantations is much smaller. 
A slit was placed in front of the NC in order to implant 
the c-Si strip without damaging the SiN^^ membrane on 
each side. The implanted area being 0.05x0.55 cm^, 55 
to 270 billion ions were implanted in each experiment. 
This corresponds to 0.1 ions per square nanometer. At 
such low fluence, there should be a very small proportion 
of collision cascade overlapping, so the experiments are 
comparable to simulations of single ion implantations. A 
delay of the order of 30 s occurred between the end of the 
implantation and the NC temperature scan, which lasts 
10 ms. Heat release is measured during the first tem- 
perature scan. Subsequent scans are used as a baseline 
experiments to account for the fact that the implanted 
and reference NCs are not identical. [201 121] 

Figure [l] shows as solid lines the heat release per unit 
temperature, dQ/dT, as a function of temperature, for 
different implantations. The amount of heat released is 
divided by the number of implanted ions. In order to 
compare with simulations below, results are scaled to 3 
keV, i.e., they are multiplied by a ratio of 3 keV divided 
by the implantation energy. The data collected for 10 
keV, 0.02 Si/nm^ implants featured a raw amplitude of 
less than 7 nJ/K above 500 K, a region where the signal 
becomes dominated by the thermal losses, therefore fea- 
turing significant noise, so a smoothed curve is presented. 
For implantations at LN, there is a rapid increase near 
200 K, followed by a slow, featureless decrease spanning 
several hundred degrees. Heat releases at similar ampli- 
tudes are observed for RT implants, but starting above 
400 K. At higher fluences, the signal becomes flat (not 
shown) instead of slowly decreasing ^ , a behaviour ob- 
served in poly-Si [23 and a-Si [5]. 

These measurements suggest that the disorder consists 
of structures complex enough that no single activation 
energy can be associated to their annealing. Indeed simi- 
lar experiments involving H implantation, where discrete 
processes occur, do show distinct peaks, associated to 
vacancy mobility, on top of a broad background signal, 
related to more complex processes [25|. 

We first compare experiments to the IV-pair model 
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FIG. 1: Heat release as a function of temperature 

starting from (a) liquid nitrogen (LN) and (b) at room 

temperature. Solid lines: Experimenta results; the NC 

signals are divided by the number of implanted ions and 

scaled to 3 keV (see text). Dashed red lines: model 

derived from k-ART after 3-keV Si ion implantation. 

Dotted black curve: simulation based on the IV-pair 

model at LN temperature. 



proposed by Marques [13 using as input the list of de- 
fects predicted by SRIM with an 8w keV implantation [9j . 
Defects are then evolved using lattice-based kinetic MC 
annealing during and after implantation with an activa- 
tion energy £^act that depends on the number n of nearest 
neighbours in the "IV-pair" state [13 . The simulation in- 
volves a 30 s waiting period at the lowest temperature to 
account for the delay between the end of the implanta- 
tion and the beginning of the NC temperature scan. A 
typical result of our implementation of this model, with a 
displacement energy of 5 eV and at liquid nitrogen tem- 
perature, is presented in Fig. [l] (a) (dotted curve) with 
the signal scaled by 3/80 to compare to experiments and 
other models. The simulation shows an isolated peak 
at 230 K, associated with the annealing of isolated IV 
pairs (n = 0) associated with activation energy of 0.43 
eV, followed by a series of peaks corresponding to the 
other activation energies (n > 0), and fails to reproduce 
experiments. Clearly, important relaxation mechanisms 
are missing from this lattice-based KMC model. 

To attempt to better capture the complex relaxation 
between the various structures forming the disorder, we 
turn to the kinetic Activation-Relaxation Technique, k- 
ART, an off-lattice self-learning kinetic Monte Carlo 
method. K-ART couples ART nouveau for the event 
search with a topological analysis tool for the catalog 
building and with a KMC algorithm for the time prop- 



agation [181 E] • By performing an extensive search for 
saddle points and fully relaxing the relevant energy barri- 
ers before each event, k-ART is able to take into account 
short and long-range elastic effects. The topological anal- 
ysis, performed with NAUTY [26 , allows a stable and 
reliable management of events even for disordered and 
complex configurations such as vacancy diffusion in Fe 
[27] and relaxation in amorphous silicon [T9|. 

We first inject a 3 keV Si atom into a 300 K 100 000- 
atom slab of Stillinger- Weber silicon with a Langevin 
bath as boundary condition perpendicular to the trajec- 
tory of the implanted atom, simulating the effect of a low- 
energy low-fluence ion implantation experiment. The cell 
is then relaxed using MD for one to ten ns. A 27 000- 
atom subpart of the cell containing the defects is then cut 
out and placed into a box with periodic boundary con- 
ditions along all direction, to eliminate surface effects. 
We generated three independent samples following this 
procedure and then launched several 300 K k-ART sim- 
ulations on each MD-produced cascade, simulating up to 
timescales of 1 s or more. 

Energy evolution for the k-ART trajectories is shown 
in Fig. |2] (a). Each relaxation simulation shows a loga- 
rithmic time dependance over many orders of magnitude 
with a slope determined by the initial disordered con- 
figuration; independent k-ART simulations starting from 
the same initial model follow a similar relaxation path- 
way. Such logarithmic relation is similar to that observed 
in heat release and relaxation experiment in other com- 
plex systems such as polymer glasses [4] . After averaging 
over the various simulations, in order to reproduce exper- 
imental measures over a large number of cascades, we find 
a slightly curved but still logarithmic overall relaxation 
behavior for the system that is robust to the addition or 
subtraction of other runs (Fig. [2] (b)). 

Before examining in more details the atomistic mecha- 
nisms leading to the macroscopic relaxation, we need to 
confirm that the simulation results reproduce the experi- 
mental data. For this, it is necessary to convert our fixed 
temperature results into constant-heating energy release. 
We consider an initial density of processes n{E, t = 0) 
that can be activated by first order kinetics. The num- 
ber of activated events with a barrier between E and 
E -h dE during the time interval between t and t + (it is 



dn{E,t)dEdt = -n{E,t)ue-^^^^ dEdt. 



(1) 



This equation can be solved at fixed (similar to the simu- 
lation) or increasing temperature (such as in NC scans). 
Following the analysis of the accepted events, and in 
agreement with previous work [6], we assume that acti- 
vated processes release a fixed amount of energy h^^ inde- 
pendent of barrier height. A direct analysis of the micro- 
scopic events indicates that the effective n(£', 0), which 
is kinetically determined, goes like n{E^ 0) ex E^ with 
—2.0 < /i < —1.2. To improve accuracy, we fit Eq. [l]to 
the averaged fixed temperature energy relaxation shown 
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FIG. 2: (a) Energy evolution for k-ART simulations. 

The lines correspond to independent simulations 

starting from three different implentation runs, (b) Red 

solid line: Average of the curves in panel (a). Green 

dashed curve: Solution of Eq. [T]at T = 300 K. Insets: 

typical configurations of defects in a 27000-atoms 
system after a 10 ns MD run (upper-right panel; there 

are 184 point defects) and after a 0.1 s k-ART 
simulation (lower- left panel; there are 94 point defects). 
Interstitials are colored in beige and vacancies in blue. 



in Fig. [2] (b) and find that the event density per interval 
log(t) decreases with increasing time, in agreement with 
the direct estimation, and that hQn{t = 0) = b?>E~^-^ . 
With this distribution, we compute the heat released 
under experimental conditions, including an initial an- 
nealing period of 30 s at the implantation temperature 
to reflect experimental conditions. Results are shown in 
Figs. [1] (a) and[l](b). The agreement between simulation- 
derived energy release and experiments both at liquid 
nitrogen (LN) and room temperature (RT) is excellent, 
especially considering that the simulations are performed 
with an empirical potential. 

The correspondence between simulations and exper- 



iments is non-trivial. A uniform density of process, 
n{E, ^ = 0), for example would lead to a flat heat release 
such as that observed in higher implantation fiuences [25] 
and a-Si [6l |23]. The origin of the exponent in the den- 
sity of processes is therefore associated with a decrease in 
the number of available barriers as defects anneal, a be- 
haviour that does not take place when amorphous pock- 
ets remain present. This decrease in complexity is ob- 
served in the microscopic evolution. After 1 ns following 
the implantation, before the k-ART simulation, the vari- 
ous models contain between 100 and 150 point defects (I 
or V) assembled into 20 to 30 clusters. Although these 
clusters vary in size from 1 to 30 point defects, most of 
them contain between 2 and 5 defects and are therefore 
better classified as defect complexes rather than amor- 
phous pockets. After 1 ms to 1 s, the various simulations 
lead to a range of configurations that consist, typically, 
of small defect complexes, with 3 to 5 point defects, and 
only a few clusters. This is illustrated in the insets of 
Fig. |2](b), that shows typical configurations after 10 ns 
and 0.1 s. 

Structural information is not sufficient to explain the 
logarithmic relaxation. Figure IsFa) shows the energy bar- 
rier separating each state selected by k-ART, aggregated 
over the simulations that reached at least 10 /is, as a 
function of time. Two features stand out. First, a quasi- 
continuous distribution of activated events is accessible in 
all time frames. Even after relaxing the system by several 
tens of eV, low-barrier events are still present and exe- 
cuted. Second, the maximum energy barrier for executed 
events increases logarithmically with time. 

The quasi-continuous barrier distribution suggests that 
the system can be kinetically-limited by configurational 
entropy. Most of the low-energy barriers connect fiicker- 
ing states that do not lead to structural evolution. The 
logarithmically growing maximum energy barriers, for 
their part, indicate that the structural evolution is also 
energy-limited and it is this interplay that generates a 
logarithmic energy decay. 

The link between energy relaxation and kinetics is 
given by the high-energy barriers. The inset of Fig. [31(a) 
shows the heat released by the events with an energy bar- 
rier above the 90*^ percentile, i.e. those above the straight 
line. These do not generally lead to large drops in poten- 
tial energy and many even lead to a higher-energy state, 
in agreement with recent observations that the forward 
and reverse energy barriers are totally uncorrelated [6 ; 
the energy barrier of events that lead directly to a sig- 
nificant relaxation are distributed evenly throughout all 
the selected events (Fig. |3](b)). Our simulations indicate 
that the annealing of the implanted c-Si is not system- 
atically kinetically limited by the relaxation events, but 
limited by events that allow the system to leave a region 
of the configurational space where the potential energy 
landscape is essentially fiat to reach another region where 
a large drop in potential energy is accessible. It can be 
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FIG. 3: (a) Energy barrier of all executed k-ART 

events. Points above the black line correspond to the 10 

% highest energy barriers in each time-frame. Inset: the 

potential energy asymmetry (final energy minus initial 

energy) of the executed events with an energy barrier 

lying in the top 10 % of the barriers plotted in Fig. [s] 

(a) for each time- frame, (b) Energy barrier of the 

k-ART events that release at least 0.5 eV of heat. 



described as a two-step replenish and relax process that 
explains the explains the logarithmic relaxation. High 
energy barrier events, which are mostly reconfiguration 
events, do not directly lead to a low-energy structure 
but rather unlock the system, open new low-energy path- 
ways and replenish the basin of available energy-releasing 
events. These events are associated with an almost con- 
tinuous and time-independent energy barrier distribution 
that refiects the complexity of the defects themselves and 
the impact of long-range elastic deformations. 
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